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We study the properties of surface solitons generated at the edge of a semi-infinite photonic lattice in nonlinear 
quadratic media, namely two-color surface lattice solitons. We analyze the impact of phase mismatch on 
existence and stability of surface modes, and find novel classes of two-color twisted surface solitons which are 
stable in a large domain of their existence. © 2008 Optical Society of America 
OCIS codes: 190.0190, 190.6135 
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Surface modes appear as a special type of waves lo- 
calized near an interface separating two different me- 
dia. In optics, linear electromagnetic surface waves are 
known to exist at an interface separating homogeneous 
and periodic dielectric media [1]. Recently, the interest 
in the study of electromagnetic surface waves has been 
renewed, and it was shown theoretically [2-4] and exper- 
imentally [5-8] that nonlincarity-induced self-trapping 
of light may become possible near the edge of a one- 
dimensional waveguide array leading to the formation 
of discrete surface solitons (see also the review pa- 
pers [9, 10]). In particular, it was found that the self- 
trapped surface modes acquire some novel properties 
different from those of the discrete solitons in infinite 
lattices: discrete surface states can only exist above a 
certain threshold power and, for the same value of the 
power, up to two different surface modes can exist si- 
multaneously. This can be understood as discrete optical 
solitons [11] localized near the surface and the action of 
a repulsive force from the surface [3] . 

Surface solitons are usually studied for cubic or sat- 
urable nonlinear media. However, Siviloglou et al. [6] re- 
ported on the first observation of discrete quadratic sur- 
face solitons in periodically poled lithium niobate waveg- 
uide arrays. By operating on either side of the phase- 
matching condition and using the cascading nonlincar- 
ity, they observed both in-phase and staggered discrete 
surface solitons. They also employed a discrete model 
with decoupled waveguides at the second harmonics to 
model some of the effects observed experimentally 

The purpose of this Letter is two- fold. First, we extend 
substantially the earlier theoretical analysis performed 
by Siviloglou et al. [6] and study, for the first time to 
our knowledge, two-color quadratic surface solitons in a 
continuum model with a truncated periodic potential. 
We analyze the effect of the mismatch on the existence, 
stability, and generation of surface solitons located at 
the edge of a semi-infinite waveguide array in nonlinear 
quadratic media. Second, we reveal the existence of novel 
classes of surface solitons which are stable in a large do- 
main of their existence. 

We consider propagation of light at the interface of 




Fig. 1. (Color online) Profiles of two-color odd surface 
solitons with (a) b\ = 1.5 and (b) b\ — 2. Profiles of 
two-color twisted surface solitons with (c) b\ = 1.75 and 
(d) b\ = 2.2. Black and red curves show the profiles of 
FF and SH fields, respectively. Lattice depth p = 1, and 
phase mismatch j3 = 0. In the white regions R(x) < 1, 
while in the gray regions R(x) > 1. 

semi-infinite lattice imprinted in quadratic nonlinear me- 
dia, which involves the interaction between fundamental 
frequency (FF) and second-harmonic (SH) waves. Light 
propagation is described by the following coupled non- 
linear equations [12,13] 



dz 



<h d 2 qi 
2 dx 2 
d 2 d 2 q 2 
2 dx 2 



qlq 2 exp(-i(3z) - pR(x)qi 
q 2 exp(if3z) - 2pR(x)q 2 . (1) 



where q\ and q 2 represent the normalized complex ampli- 
tudes of the FF and SH fields, x and z stand for the nor- 
malized transverse and longitudinal coordinates, respec- 
tively, P is the phase mismatch, and d± = — 1, d 2 = —0.5; 
p is the lattice depth; the function R(x) = at x < 
and R{x) = 1 — cos(Kx) at x > describes the profile 
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Fig. 2. (Color online)(a) Power of surface solitons vs. 
propagation constant for different phase mismatches, (b) 
Cutoff of the propagation constant b\ for surface solitons 
as a function of depth at (3 = 0, Inset shows the cutoff vs. 
phase mismatch at p — I. (c) Power of twisted surface 
solitons vs. bi for different phase mismatches, (d) Cutoff 
of the propagation constant for twisted surface solitons 
as a function of lattice depth at (3 = 0. Inset shows the 
cutoff vs. the phase mismatch at p = 1. 

of a truncated periodic lattice with modulation K . The 
system of Eq. ([1]) admits several conserved quantities in- 
cluding the power P = /^(kil 2 + \q 2 \ 2 )dx. 

The stationary solutions for the lattice-supported sur- 
face solitons can be found in the form q\^ 2 (x,z) = 
ui, 2 {x)exp(ibi^ 2 z), where u\^ 2 (x) are real functions, and 
&i,2 arc real propagation constants satisfying 62 = P + 
2&i. Families of surface solitons arc determined by the 
propagation constant b\, the lattice depth p, and the 
phase mismatch (3. For simplicity, we set the modula- 
tion parameter K = 4. To analyze stability we examine 
perturbed solutions qi i2 (x,z) = [ui >2 {x) + Ui. 2 (x,z) + 
iVi,2{x, z)]exp(ibiz), where real parts Ui t2 and imagi- 
nary parts V1.2 of perturbation can grow with complex 
rate S. Linearization of Eq. (JT|) around Ui t2 yields the 
eigenvalue problem 



5U X 



di d 2 V x 



{u{V 2 - U2V1) - pRVi + biVi 



2 dx 2 
6V i = -%^TT + {uiU 2 + u 2 U 1 )+pRU 1 - hUi 



SU 2 



SV 2 



2 dx 2 
d 2 d 2 V 2 
2 dx 2 
d 2 d 2 U 2 
2 dx 2 



2«iVi - 2pRV 2 + b 2 V 2 



2u 1 U 1 + 2pRU 2 -b 2 U 2 . (2) 



which we solve numerically and find the growth rate S. 
Figures [lja,b) show the examples of the simplest two- 
color surface lattice solitons, the so-called odd solitons. 
Such surface modes reside in the first lattice channel, 
where both FF and SH fields reach their maxima. Due 





mismatch 

Fig. 3. (a) Instability domain (shaded) for twisted sur- 
face solitons on the (&i , (3) plane at p = 1. (b) Real part 
of the instability growth rate for twisted surface solitons 
vs. propagation constant at (3 = and p = 0.5. 



to the lattice truncation, the solitons have asymmetric 
profiles at lower power [Fig. [TJa)]. There exists a lower 
cutoff (&co) of the propagation constant for the existence 
of odd solitons. The power of odd solitons is a nonmono- 
tonic function of the propagation constant, and there is a 
narrow region close to the cutoff where the power depen- 
dence changes its slope, dP/dbi < [Fig. Ufa)]. We find 
that the cutoff 6co is a monotonically increasing func- 
tion of lattice depth p, and it is a decreasing function 
of the phase mismatch [Fig. [2jb)] . It should be noted 
that the critical power for the existence of two-color sur- 
face lattice solitons depends on the phase mismatch, it 
reaches the minimum value for the exact phase matching 
(/? = 0). Linear stability analysis reveals that these odd 
surface solitons are stable almost in the whole domain 
of their existence except a very narrow region near the 
cutoff b co where the exponential instability develops. Di- 
rect numerical simulations of the model |T]) confirm the 
results of the linear stability analysis. 

In addition to the simplest solitons described above, 
we find various families of higher-order surface lattice 
solitons, which can be viewed as the combination of sev- 
eral in-phase and out-of-phase odd solitons. Here we only 
focus on the case where the FF field features the out- 
of phase combinations because the modes with in-phase 
combinations arc unstable. Being reminiscent to its dis- 
crete counterparts such as twisted localized modes, here 
we term the mode with out-of phase combination as 
twisted surface lattice soliton. Typical profiles of twisted 
modes residing at the edge of a semi-infinite lattice are 
shown in Figs. [Uc,d). Similar to the properties of odd 
solitons, the power of twisted surface solitons is a non- 
monotonic function of the propagation constant, and 
there exists a narrow region close to the cutoff where 
dP/dbi < [Fig. (He)]. However, one should note that 
the cutoff b C o is a nonmonotonic function of the lat- 
tice depth p, as shown in Fig.[2jd). Likewise, the critical 
power for the existence of twisted solitons reaches the 
minimum value at (3 = . 

Importantly, we find that the twisted surface lattice 
solitons become completely stable when the power ex- 
ceeds a certain threshold value [Fig. [3^a)], while the in- 
stability domain expands with the growth of the absolute 
value of phase mismatch (3. Figure [3]Jb) shows the real 




Fig. 4. (Color online) Profiles of two-color surface soli- 
tons in the second channel: (a) odd one with 61 = 1.5 
and (b) twisted one with b\ — 1.75. The black and red 
curves show the profiles of FF and SH fields, respectively. 
Power of (c) odd and (d) twisted surface solitons in the 
first (blue) and second (black) channel of the lattices. 
Here lattice depth p = 1, and phase matching (3 = 0. 



part of the perturbation growth rate versus the propa- 
gation constant, where the left part corresponds to the 
exponential instability, while in the right part solitons 
suffer from oscillatory instabilities. Results from linear 
stability analysis are confirmed by the direct numerical 
simulation of the model ((T|) . 

We also study the properties of two-color surface 
modes which are shifted from the lattice edge, similar 
to the analysis performed earlier for the discrete cubic 
model [3]. Figurcs0Ja,b) show some illustrative examples 
of odd and twisted states residing in the second channel 
of the semi-infinite waveguides from the interface. Such 
modes require much lower critical power comparing with 
surface modes residing at the interface [Fig. Hfc,d)]. Lin- 
ear stability analysis reveals that two-color surface soli- 
tons in the second channel are stable when their power 
exceeds a certain threshold. 

To address the issue of excitation of two-color surface 
solitons, we perform a comprehensive study of the soliton 
generation by two Gaussian beams qi t 2 = ai,2exp(— x 2 ). 
As shown in Figs. Efa-d), for lower powers the beams 
just experience repulsion from the surface and diffraction 
[(a)]. For high enough powers, surface solitons can be 
excited either from only FF field [(b)] or both FF and 
SH fields [(c)] with a proper phase matching. Two-color 
surface solitons can be formed for other conditions (/3 = 
—3), but for higher input powers [(d)]. 

In conclusion, we have analyzed the existence, stabil- 
ity, and generation of two-color quadratic surface solitons 
in a continuum model with a truncated periodic poten- 
tial, and also revealed the existence of novel classes of 
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Fig. 5. Excitation of two-color surface solitons from 
Gaussian input beams: (a) a\ = 1.5, a-i = 0; (b) a\ = 2.8, 
<22 = 0; (c) a\ = 1.5, 02 = 1; and (d) a\ — 2, a-i = 2.5. 
Here (a-c) /3 = 0, and (d) (3 = —3. In all cases, p = 1 
and only the SH field is shown. 



stable parametrically coupled surface states. 
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